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Abstract. The average avalanche size can be calculated exactly in a number of models of self-organised 
criticality (SOC). While the calculation is straight-forward in one dimension, it is more involved in higher 
dimensions and further complicated by the presence of different boundary conditions and different forms 
of external driving. Amplitudes of the leading order are determined analytically and evaluated to obtain 
analytical references for numerical work. A subtle link exists between the procedure to calculate the average 
avalanche size and the field theory of SOC. 
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1 Introduction 

The average avalanche size in models of self-organised crit- 
icality (SOC) [1] is one of the few observables whose scal- 
ing frequently is known exactly even in non-trivial cases. 
In numerical simulations, it often plays the role of a bench- 
mark for convergence [2,3]. In the one-dimensional ver- 
sion of many models, even the amplitude of the average 
avalanche size is known exactly, as well as the corrections 
to scaling [4]. In the following, exact results for the av- 
erage avalanche size are collated and extended to higher 
dimensions and more complicated boundary conditions. 

In a number of conservative SOC models [1,5], parti- 
cles (height units or slope units) perform a random walk 
from the point of being added to the system until they 
leave it. In the Manna model [6] particles move indepen- 
dently from site to site, so that their trajectory is exactly a 
random walk. In that model, particles are added at (nor- 
mally randomly and independently chosen) sites by the 
so-called external drive. If the number of particles at a 
site exceeds 1, all particles are being redistributed inde- 
pendently and randomly to the nearest neighbours. This 
process is repeated until no height exceeds the threshold of 
1 anymore. Sites that do (and the particles located there) 
are called active. The totality of redistributions or top- 
plings constitutes an avalanche. A complete separation of 
the time scales of driving and relaxation is achieved by 
driving only when no site (or particle) is active. A parti- 
cle may rest for very long times until being moved again, 
but while it is moving, it performs a random walk in a 
time that advances only as long as the particle is active 
(conditional time scale). 



Even in systems where particle movements are not in- 
dependent, such as the BTW [1] and the Oslo Models 
[7], where particles are redistributed evenly among near- 
est neighbours, it has been noted that particles follow 
random-walker trajectories, because the entire ensemble 
of possible paths is being generated as sites topple. 

It has therefore been noted several times that the aver- 
age avalanche size in the Manna Model is essentially given 
by the average escape time of a random walker, e.g. [8,9]. 
While particles in the Manna Model describe trajectories 
of a random walker, each of their moves from one site to 
a neighbouring site is caused by a toppling. In fact, in 
the Abelian version [10] considered in the following, each 
toppling causes two particles (2d particles in the BTW 
and Oslo Models on hypercubic lattices 1 ) to move and so 
the average number of topplings per particle added, which 
is exactly the average avalanche size, is equal to half the 
average number of moves each particle makes until its de- 
parture from the system. 

The number of charges a particle causes during its life- 
time (i.e. the number of times a particle arrives at a site 
until it leaves the system) , is exactly equal to the number 
of moves it makes; while the initial deposition represents 
a charge, but not a move, the final move (off the system) 
does not cause a charge. 

As opposed to higher moments, the average avalanche 
size can be calculated because it does not require any in- 



More generally, in BTW and Oslo Models, q particles top- 
ple, where q is the coordination number of the lattice in the 
bulk. 
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formation about the collective toppling of particles. 2 It 
is merely a matter of stationarity and conservation. The 
former is important because only at stationarity the av- 
erage avalanche size can be determined as the number of 
topplings per particle exiting by averaging over so many 
avalanches that the vast majority of particles added have 
left the system. Conservation is important for two rea- 
sons. Firstly, particles should not disappear by interaction, 
which cannot be accounted for in this simple approach. 
Secondly, each and every toppling must count towards an 
avalanche. 

In the following, the average avalanche size is calcu- 
lated for hypercubic systems in arbitrary dimensions (but 
see Section 5). First, it is calculated for a one-dimensional 
"lattice" with two open boundaries. The result is then 
generalised to the scaling in arbitrary dimensions. Doing 
this exactly and on the lattice is a difficult undertaking [9] , 
but the aim of the following is to determine the leading 
order amplitudes. 3 After taking the continuum limit, they 
are calculated for a variety of boundary conditions. Some 
special cases are discussed. Finally, the result is related to 
some recent field theoretic insights. 



2 One dimension 

In one dimension, the average number of moves can be cal- 
culated fairly easily for a variety of boundary conditions. 
For brevity, I focus on two open boundaries (i.e. particles 
leave the system if a toppling site attempts to deposit a 
particle on an "outside" site). If x is the site a particle is 
added to by the external drive, then the average number 
of moves m(xo;L) the particle makes until its departure 
is given by [11, 12] 



m(x; L) = 1 + 



m(x + 1;L) + m(x — 1; L) 



(1) 



where the open boundaries are implemented by imposing 
m(0;L) = m(L + 1;L) = 0, i.e. a Dirichlct boundary 
condition. Rearranging terms produces a Poisson equation 
on the lattice, whose solution is a simple quadratic, 



m(x; L) = x(L + 1 — x) . 



(2) 



Summing over the uniform drive (i.e. xq uniformly and 
randomly taken from {1, . . . , L}) gives 



m(L) = lX>( I; L)= <^'f + 2 > 



(3) 



X = l 



and thus the expectation of the avalanche size (first mo- 
ment) is exactly [13] 

M-l-W-™.!'. (4) 



2 In contrast, the present approach does not allow the cal- 
culation of the average avalanche size in the ensemble of 
avalanches with non- vanishing size. 

3 In the following, when quoting results to leading order the 
equality sign ~ will be used. 



2.1 Generalisations 

In higher dimensions, the scaling (s) oc L 2 persists, which 
is of course just the usual escape time of a random walker: 
It explores the distance L within L 2 moves. This argu- 
ment can be made more rigorous by noting that if the 
survival probability after t moves (i.e. the probability of 
the random walker not having reached an open boundary) 
is o~(t, L) in one dimension (for the sake of simplicity, this 
is the probability averaged over the uniform drive), then 
in higher dimensions d that probability is simply a(t, L) d , 
because of the independence of the d directions of possi- 
ble displacement and the hypercubic nature of the bound- 
aries. 4 In the continuum limit, t is better interpreted as 
a time, rather than a number of enforced moves. The av- 
erage residence time in d dimension, equal to the average 
time to escape £<j, is thus 5 

poo 

2{s) d (L)~e d (L) = / dta(t,L) d 
Jo 

= f< < 5 > 

where — -^u(t, L) d is the probability density of escaping at 
time t. Its structure reflects the fact that the movement in 
the d spatial directions is independent; — ^cr(i, L) is the 
probability density to escape at time t in one direction, of 
which there are d (choices), and a(t 1 L) d ~ 1 is the proba- 
bility to stay within bounds in the remaining d— 1 dimen- 
sions. Here and in the following, the factor 2 in front of (s) 
(on the left of Eq. (5)) is retained, acting as a reminder 
of its origin as the number of particles redistributed in 
each toppling (the avalanche size s being measured by the 
number of topplings). In the BTW and the Oslo Models, 
that factor 2 has to be replaced by the coordination num- 
ber of the lattice, 2d for a hypercubic one with nearest 
neighbour interaction. 

Because cr(t,L) is, by dimensional consistency, bound 
to be the dimensionless function a(t/L 2 , 1) it follows that 
Ed oc L 2 , in line with the view that the trajectory of a ran- 
dom walker is essentially a two-dimensional object [14]. 
Claiming that t/L 2 is dimensionless means being some- 
what cavalier about the dimension of the diffusion con- 
stant D, which in the present context relates time and 
number of moves. If the walker takes, in each time step, 
one step in any of the d spatial directions, the variance of 
its displacement is 1. The diffusion constant, on the other 
hand, is half the variance of the displacement in each (in- 
dependent) spatial direction per time, so that 2Dd = 1 on 
hypercubic lattices. There is thus a slight conceptual dif- 
ference between the active particles in the Manna Model 
on the one hand, which are forced to move to one of their 
nearest neighbours, and a random walker with a certain 



4 If the boundaries are shaped or structured then the survival 
in one direction depends on the coordinate in the other. Results 
for that case can be found in [9]. 

5 The ~ sign applies as Ed(L)/2 is a continuum approxima- 
tion of (s} d (L), yet ea(L), in the continuum, itself is calculated 
exactly. 
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diffusion constant on the other, which is subject to ran- 
dom motion in each spatial direction independently. 

The survival probability can be calculated quite eas- 
ily, noting that the normalised eigenfunctions of d 2 
with Dirichlet boundary conditions in one dimension are 

y/2/ L sm(xq n ) with q n — rnr/L, where n = 1, 2, With 

periodic boundary conditions, they are exp(xq n ) with 
q n = 2nir/L and any integer n, including and negative 
integers, n £ Z. As it will turn out below, given the self- 
adjoint operator d 2 , it is the presence or absence of the 
zero mode, i.e. the constant eigenfunction with eigenvalue 
0, which decides over conservation or dissipation and the 
structure of the resulting equation for (s). 

In one dimension, the probability density function 
(PDF) of a particle under Brownian Motion started at 
x with diffusion constant D on an interval with open 
boundaries at and L is thus 



V {x, t; x , L) = — ^ sin(xq n ) sm(x q n )e Dq " 



(6) 



n=l 



Since the motion in the different directions is independent, 
the PDF in higher dimensions is a product of Eq. (6). The 
expected residence time is given by the integral over time 
and space, in one dimension 



pL />oo 2 00 

dx dt V (x, t; x , L) = — 

Jo Jo L „~r^ 



n— l.odd 



2 ■ f A 1 
Qn Dqi 



(7) 

where the constraint of n = 1, 3, 5, . . . in the sum being 
odd comes from the integral of ain.(xq n ), which gives 2/q n 
for odd n and otherwise. For uniform drive the escape 
time is given by the integral j^dxa {l/L) of Eq. (7). The 
survival probability, a(t,L), on the other hand, is given 

by 



a(t,L) = — ( dx f dxV (x, t; x , L) 
L Jo Jo 



L 2 V ql [ ' 



and therefore 



2 (s) d (L) ~ e d (L) 



_2_y ^ ±j_ 

7 n,m,.. -l.odd Hn y ™ 



"£>(<£ + «& + ...) 



(9) 



where the sum runs over d different indeces. Dhar's re- 
sult for the lattice in d — 2 [15, Eq. 21] is recovered by 
approximating cot(7m/(2L + 1)) w 1/ sin(7rn/(2L + 1)) w 
(2L + l)/(7rn) for large L. By comparison with his results 
it is clear that in general, on hypercubic lattices the con- 
fluent singularities in the finite size scaling of (s) are L 1 , 
L° etc. 



2.2 One dimension again 

In the following, a few particular results deriving from 
Eq. (9) are highlighted. In one dimension, 



, . / r \ 2 V- 4L 2 2L 2 



n— l,odd 



t-2 

T < 10 > 



using E~=i,odd V" 4 = tt 4 /96 [16, Sees. 1.471 and 1.647], 
consistent with Eq. (4). Sums of this type frequently oc- 
cur in finite temperature field theory under the label of 
Matsubara sums [17]. The latter is associated with the 
technique of representing the sum as one over residues, 



oo 1 1 
2 V — = — 



1 



2m J c z 4 1 + exp {ittz) 



(11) 



where the contour C (see Figure 1) encircles each (simple) 
pole of —tn/(l + exp(i7T2;)), which are located at z = n 
and z = —n (n odd; the parity symmetry is the origin 
of the factor 2 on the left) and have residue 1. Merging 
the contours for z = q n and z = —q n and deforming the 
resulting two contours to enclose the single pole of order 
4 at z — produces the desired result, as the contour has 
negative orientation and the residue is — 7r 4 /48. 

It is instructive to attempt to recover Eq. (2), which is 
twice the avalanche size for a system driven at site x = a. 
In that case, the uniform drive, J dxo (l/L), has to be 

replaced by a single source at a, i.e. J^dxQd(a — xq), so 
that the average avalanche size for a system driven at a is 



2 <*>i,„ (£) 




2L fmra\ 2L 2 
— sin — — -3-2 



sin 



71=1, odd 



/ mm \ 1 



(12) 



Clearly the terms in the sum contribute significantly less 
for large n. For small n and large L, the sin may be ap- 
proximated by its argument, producing 



L* v-^ • /nna\ 1 1 ^ 



a 

n=l,odd 

and thus 



(13) 



n— l.odd 



2(s) ha (L)~La (14) 

as in Eq. (2). One may be tempted to consider the sum as 
Ricmann sum with mesh l/L, 



E 

n=l,odd 



sin 



/■nna\ 1 



du sin(u) 



1 



with dummy variable u = anir/L. In the last step, the 
integrand has been approximated by 1/u 2 valid at small 
u. In this approximation, the avalanche size is 

2 <*>i,„ (£)«-^r ' ( 16 ) 
a rather poor approximation compared to Eq. (14). 
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(a) Initial arrangement of simple poles to evaluate sum 
Eq. (11). 



X 



X) « (X 



X 



(b) Deformation of the initial contour to form two big ones. 



X 



X 



(c) Deformation of the contour enclosing a different pole. 

Fig. 1. Arrangement of the poles in Eq. (11) in the complex 
plane and subsequent deformation of the contour, (a) The sum 
Eq. (11) can be performed in a contour integral, by evaluating 
1/z 4 when calculating the residues at the poles of a suitable 
function, indicated by crosses. All circles (with arrows indicat- 
ing the direction) together make up the (initial) contour C in 
Eq. (11). (b) The contour can be deformed by merging the in- 
dividual paths. The pole of the factor 1/z 4 (filled circle) has to 
be avoided. The dotted lines indicate the form of the contour 
for z — > ±oo. (c) Because the integrand drops off sufficiently 
fast in large arguments, the contours can be joined up (dotted 
line) differently, thereby generating a path around the pole of 
1/z 4 with opposite orientation. 



numerical estimate Eq. (18), the average avalanche size in 
two dimensions with uniform driving in the bulk and open 
boundaries is 



2( S ) 2 (L)^C 2 



64 



(19) 



with C 2 = 0.5279266525115576573... and D = 1/4, so 
that 



(s) 2 (L) ~ 0.070288507477576867. . .L 2 



(20) 



Extensions of the result above to systems with non- 
unity aspect ratio r = L x /L y are straight forward. The 
sum to be performed is 



2 (s) 2 (L x , L y ) — j^-qL x L v 
and thus 



1 



^—f n 2 m 2 (rn 2 +r 1 m 2 ) 

n,m— l,odd v ! 

(21) 



2 (s) 2 (L x ,L y ) 
_ 64 

- Jj^E x v 



768r 



E 

771=1, odd 



7rtanh(7T77ir/2) 
4r 2 m 5 



(22) 



which according to Eq. (21) is invariant under a change of 
r to r~ l . While this is not at all obvious in Eq. (22), the 
sum is somewhat reminiscent of that in [16, Sec 1.471.3]. In 
the limit of large r, the first term in the curly bracket dom- 
inates, producing 2 (s) 2 (L x ,L y ) = L 2 /(12D), as the sys- 
tem is essentially one-dimensional (except for D=l/(2d), 
due to the additional degree of freedom). For very small r 
the sum may be treated as a Riemann sum. 



3 Two dimensions 

In two dimensions, the same procedures can be followed. 
For uniform driving, the key sum to perform is 



E 



i 



n,m— l,odd 



n 2 m 2 (n 2 + m 2 ) 

E7T 2 7rtanh(-7rm/2) 
8m 1 4^5 

m— l,odd 

7r 6 ^ 7rtanh(7T77j/2) 

~ 768 ~ ^ 



m— l,odd 



Am 5 



(17) 



While the author was unable to determine the last sum 
(relevant literature [18,16,19]), it converges extremely 
quickly because of the large power of m in the denomi- 
nator and because tanh very quickly approaches 1. With 
the help of Mathematica [20] , one finds 



E 

m— l.odd 



tanh(7rm/2) 



0.9216754342259668... (18) 



to be compared to tanh(7r/2) = 0.91715233... and 
tanh(7r/2) + tanh(37r/2)/3 5 = 0.92126689.... Using the 



3.1 Mixed boundary conditions 

If any of the boundaries is periodically closed or made 
reflecting, the dimension (i.e. this degree of freedom) ef- 
fectively disappears from the problem, i.e. the avalanche 
size is essentially that of a one-dimensional system. The 
only trace that remains of the "closed dimension" is hid- 
den in the diffusion constant, which is D = l/(2d), so 
that 

<«>2,cyl = Y +0(L) (23) 

for cylindrical boundary conditions on square lattices, d = 
2. In Eq. (23) the sub-leading terms are indicated as well, 
because Eq. (4) remains exact, i.e. 



<*>2 



(L + l)(L + 2) 



'2,cyl - 6 

and obviously in higher dimensions 

_ (L + l)(L + 2)d 



cyl 



12 



(24) 



(25) 



is the exact expected avalanche size on a hypercubic lat- 
tice if only one direction remains open, while d — 1 direc- 
tions are either periodically closed or reflecting (or, more 
generally, produce a spectrum containing 0). 
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The technical reason for the simplicity of the results 
with cylindrical boundary conditions is the presence of 
only one sum. The other sums do not occur because the 
integration over the entire system as well as the integra- 
tion over uniform drive effectively projects the eigenfunc- 
tions of the PDF on a constant, using the scalar product 
with constant weight. Under that weight V 2 is self-adjoint 
and the constant has eigenvalue 0, provided it is an eigen- 
function (which is decided by the boundary conditions). 
In that case, the <?o = mode is selected in the sum, so 
that neither any factor 1 /q n appears nor a contribution in 



.). The mode with eigenvalue does not 



decay in time, i.e. it is conserved. A boundary condition 
that leads to conservation is thus expected to possess such 
a 0-mode. 

The same type of argument therefore applies in higher 
dimensions. For example, when applying periodic bound- 
ary conditions to d — 2 directions in hypercubic lattices 
with d > 2, the average avalanche size is essentially that 
of a two-dimensional system, except for D being changed 
to£> = l/(2d). 

If individual boundaries have mixed conditions, cal- 
culations become drastically more complicated. An inter- 
esting example is a setup where all boundaries of a two- 
dimensional lattice are reflecting except for a narrow hole 
of fixed size h from where all particles added have to es- 
cape. Rather counter-intuitively, the scaling of the escape 
time in L is not very different from the scaling on the 
open lattice, even when the size of the hole is kept finite 
and fixed as L is increased. This is surprising, as the par- 
ticles need to "find the narrow exit" in an increasingly 
large system — given the presence of an additional length 
scale (the size of the hole) the scaling of the escape time is 
no longer determined by dimensional consistency. On the 
other hand, one may argue that the situation is not much 
different from a one-dimensional lattice, where the size of 
the exit remains constant as well. 

On the basis of published results on the narrow escape 
problem [21-23], one finds 



2(a) 



2, narrow 



(L,h) 



% ( ln 



O((L/h) )\ (26) 



where h is the fixed size of the hole adjacent to a corner 
and D = 1/4 is the diffusion constant. That the size of 
the hole enters only very weakly, suggests that the escape 
time is essentially determined by the time it takes for the 
particle to explore the entire lattice, rather than the size 
of the exit hole. In fact, in dimensions d > 2 a random 
walker can be thought of as exploring a convoluted two- 
dimensional surface with the area covered (number of dis- 
tinct sites visited) increasing essentially linearly in time, 6 
as if it was never revisiting a site. Strictly, this holds only 
in dimensions strictly greater than two. In two dimensions 
and less, every site is returned to infinitely often. 

As a final special case in two dimensions, I consider a 
system driven at a site with fixed coordinates while the 



The walk in one dimension can be interpreted as a projec- 
tion from two dimensions. 



system size is increased. After the considerations above, 
it is clear that fixing d < d coordinates will lead to a 
scaling (s) oc L 2 ~ d " for g? < 2 in any dimension d. The 
case c?o = 2 is special because logarithmic scaling is ex- 
pected (whereas a constant average avalanche size occurs 
for g? > 2). Again, using a Riemann sum introduces un- 
controllable errors which are exacerbated by the sensitiv- 
ity of the expected logarithm to small corrections. 

The average avalanche size in two dimensions with 
fixed driving position is given by 



^ ( s )2,fixcd 

1^2 22 
- ~jj 2^ —sm(q n x ) — jsm(q m y ) 

n,m— l,odd 



Qn L 



im ql + q 2 m 



(27) 



where sin(n7ra;o / L) can be approximated by its argument 
as the terms in the sum vanish at least like 1/n 3 in large 
n. For small n the resulting sum is divergent in the upper 
limit, which has to be replaced by the ultraviolet cutoff 
L/a with lattice spacing a, 



2(a) 



2, fixed 



TT 2 D 



L/a 

E 



i 



n 2 + in 2 



(28) 



The final result hinges on the last sum. One of the summa- 
tions can be performed beyond the upper cutoff without 
causing a divergence. The resulting summation involves a 
term of the form tanh(7rm/2)/m, which may be approxi- 
mated by 1/m and thus the sum by (w/8) In(L/(2a)), so 
that 

2<s> 2 , fix cd«^ln(V(2 fl )). (29) 

The role of the upper cutoff becomes clearer in the case 
da > 2, for example fixing the driving position on a three- 
dimensional lattice. The reason why the expected escape 
time remains finite even in the thermodynamic limit is 
because within a finite time the random walker, attempts 
to travel beyond the finite distance to one of the open 
boundaries, thus leaving the lattice. Without a finite lat- 
tice spacing, the number of "hops" to the open boundary, 
however, diverges. The difference between thermodynamic 
and continuum limit is that absolute distances correspond 
to a fixed number of hops in the former, but not in the 
latter. From a physical point of view, there is in fact no 
other difference between the two. 



4 Higher dimensions 

In higher dimensions the calculation of the relevant sums 
becomes increasingly computationally demanding. The 
expected avalanche size for homogeneous drive in a d di- 
mensional hypercubic system with open boundaries gen- 
erally is according to Eq. (9). 



(30) 
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where D = l/(2d) has been used and 

1 1 



then 



c d = J2 



,„ 2 _n,=o nti(2n, + I) 2 £Li(2^ + I) 2 



(31) 



One of the summations can always be carried out, 
Eq. (17). Keeping only the two lowest order terms in the 
resulting sum produces a recurrence relation for d > 1, 



4(rf- 1)( 3 / 2 
7r(d- 1) 



tanh 



and C\ = 7r 4 /96 exactly. Table 1 contains the numerical 
evaluation of the constants Cd according to Eq. (31) to- 
gether with the approximation Eq. (32). The amplitude in 
the last column are well consistent with recent numerical 
results on the Manna Model [2,3]. 



5 Arbitrary Adjacency 

Eq. (1) points to a more general procedure to calculate 
the expected number of moves to escape from the lattice. 
If |m) is a vector whose components rm are the expected 
escape times starting from site i and A is the weighted 
adjacency matrix (closely related to Dhar's toppling ma- 
trix [15], also discussed by Stapleton [12]), proportional 
to the lattice Laplacian, containing An = — 1 across the 
diagonal and A t j being the probability of i discharging to 
j (i.e. Aij = l/(2d) on hypercubic lattices), 7 then 



-\l)=A\m) 



(33) 



where |1) is a column of ones. Dissipation at boundary 
sites is implemented by J2j Aij < 0, while ^ . A^ = 
at (conservative) bulk sites. The presence of the non- 
conservative sites means that |1) is not an eigenvector, 
in fact A\l) is a vector with components that are for 
each conservative (bulk) site and negative for all dissipa- 
tive (boundary) sites. If (d\ is a vector whose components 
di are the probability that a particle is deposited at site i 
by the external drive, with normalisation (d\l) = 1, then 



2(s) = (d\m) = -(d\A- 1 |1) 



(34) 



provided the inverse A^ 1 of A exists. If A's eigenvectors 
(e$| and |e$) (not necessarily transposed relative to each 
other, as A may be directed, i.e. not symmetric), with 
eigenvalues Aj and (ej|ej) = span a subspace contain- 
ing (d\ and |1) respectively, so that 



(d| = (ei\ 



(35a) 



2( S ) = -J2(d\e i )K 1 (e t \l) = -J2 



(36) 



For uniform drive di = 1/N in a system with N sites 
and so N (d\ = (1| is a row of ones. In that case, if A is 
symmetric Uj = Wi/N and 



1 V- (ei\iy 



N ^ A, 



(37) 



6 Relation to field theory 

There is a subtle but very important link between the 
calculations performed above and the field theory of the 
Manna model [24]. Prima facie, it may look accidental that 
the calculations for the expectation of the escape time of 
a random walker are identical to those for the expected 
activity integral. In fact, the bare propagator for the ac- 
tivity at w = (vanishing frequency, as obtained after 
Fourier transforming the time domain) is identical to that 
of the time-dependent PDF of the random walker particle. 
However, while the former describes the spreading of ac- 
tivity on the microscopic time scale of the Abelian Manna 
Model [6, 10] subject to Poissonian updates (activated ran- 
dom walkers [25]), the latter describes the movement of a 
particle on the conditional time scale, which advances only 
when the particle is not stuck on the lattice. Only on that 
time scale, an actual random walk is performed and the 
link exists between the number of moves and the residence 
time. 

In the light of the field theory, however, it is clear that 
the particle movement on the conditional time scale is ex- 
actly identical to the spreading of activity; particles mov- 
ing are active and vice versa. The fact that the average 
avalanche size can be determined by the considerations 
presented above means that the bare propagator at lo = 
is not renormalised at any order. That does not imply 
that the bare propagator is not renormalised at all, as the 
statement merely applies to lu = 0. In fact, the time de- 
pendence of the propagator is very much expected to be 
affected by interaction, fluctuations and thus renormali- 
sation, because active particles do not move freely like a 
random walker, but interact with the particles at rest. 

The reason why no renormalisation of the propaga- 
tor at lj = takes place is the same reason that allows 
the calculation of the average avalanche size in the first 
place: Conservation of particles and stationarity; in the 
stationary state and because of conservation, on average 
exactly one particle leaves the system per particle added. 
The number of moves performed by a particle during its 
residence determines the average avalanche size. 



|l) = ^Wila) 



(35b) 
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d 


Cd (numerically) 


Cd (approximation Eq. (32)) 


(d/7T 2 )(8/ir 2 ) d C d 


1 


1.0146780 . . . 


1.0146780 . . . 


0.0833333 . . . 


2 


0.5279266 . . . 


0.5282475 . . . 


0.0702885 . . . 


3 


0.3737684 . . . 


0.3749565 . . . 


0.0605054 . . . 


4 


0.3026980 . . . 


0.3055630 . . . 


0.0529579 . . . 


5 


0.2651000(1) . . . 


0.2707675 . . . 


0.0469927(4) . . . 



Table 1. The constant Cd, Eq. (31), for dimension d = 1, 2, . . . , 5. The second column shows the numerical evaluation of the 
sum (with extended double precision, summing up to 2 ■ 1000 + 1 for d = 1, 2, 3, up to 2 • 500 + 1 for d — 4 and up to 2 • 200 + 1 
for d = 5). Unless an error is stated , the digits shown display convergence. The third column is the recursive approximation 
Eq. (32). The last column is the amplitude of the leading order L of the average avalanche size, Eq. (30). 
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